LEVY SOLUTIONS OF A RANDOMLY FORCED BURGERS EQUATION 
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Abstract. We consider Burgers equation forced by a brownian in space and white noise in 
time process dtu + ^dx{u)^ = f(x, t), with E{f{x, t)f{y, s)) = ^{\x\ + \y\ — \x — y\)5(t — s) and 
we show that there are Levy processes solutions, for which we give the evolution equation of the 
characteristic exponent. In particular we give the explicit solution in the case uo{x) = 0. 



1. Introduction 

We consider the randomly forced Burgers equation dtu + ^dxiu)"^ = f{x,t), where / satisfies 
E{f{x, t)f{y, s)) — A{x, y)5{t — s) with A{x, y) = \{\x\ + \y\ — \x — y\): / is a brownian in space/ 
white noise in time process (the brownian is considered on the whole space, and is obtained with 
two independent brownians, one on the left and one on the right). Burgers equation has originally 
been introduced as a ID model of turbulence. It is certainly too crude, but it is still a good idea 
when trying to find new approaches to the "real" 3D turbulence as given by Euler equation, to 
start by looking at Burgers. In this point of view, the randomly forced Burgers equation is a crude 
simplification on the way to the description of forced turbulence. A variant of our problem was 
considered by [5, where they considered a gaussian in space/white noise in time forcing, satisfying 
E{f{x^t)f{y,s)) = A{x — y)6{t — s) and obtained all the hierarchy of n-point densities evolution 
equations. We will proceed quite differently, basing our study on the result shown by Carraro 
and Duchon in [6], and by Bertoin in [3], that Levy processes with negative jumps are statistical 
solutions (in a sense that we will recall later) of the Burgers equation without any forcing. This is 
actually a particular case of the fact that the Burgers equation conserves the Markovian (in space) 
property of processes ([7]). We will show here in the same spirit that the brownian in space/ white 
noise in time forcing allows to keep the Levy property. Moreover, our method allows us to write 
explicitly the exponent of the Levy process that one gets when starting from uq{x) = 0. 



2. Statistical solutions 

We will closely follow [6] (see also jl3]). Let E be the space of cadlag real functions. We will 
call C{E) the smallest cr-algebra such that for each x £ M, u i— > u{x) is measurable, and C'{E) the 
smallest cr-algebra such that for each (z, y) e M x M, u u(x) — u{y) is measurable. Let Vq be 
the set of real C°° functions v with compact support such that Jj^ v{x)dx = 0. A probability /i on 
{E,C{E)) is then characterized by its characteristic function 



V eV t-^ fl{v) := / exp [i / u{x)v{x)dx] dii{u) 
J E Js. 

whereas a probability n on {E,C'E)) is characterized by its characteristic function v G Vq i-^ P-{v)- 
Let uo '■ (il,yl, P) —fEhea random process, defined on some probability space, and let 
fiQ : C{E) — !■ [0, 1] denote its probability law : uq will be our initial condition. We will let u 
evolve according to the non forced Burgers equation, obtaining thus a family of processes ut where 
ut{x) = u{x,t) is the solution of Burgers at time t. Then it is not difficult to check that the 
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characteristic function /tj of Ut verifies : 

dtfLt{v) = dt{ exp [z / u{x)v{x)dx]dnt{u) 



E 



E 



9({exp [i / u{x,t)v{x)dx\} diia{uQ) 



exp j u{x,t)v{x)dx)](dt[i j u{x,t)v{x)dx]) d^j.o{uQ) 

1 



(f) = / exp [i / u{x)v{x)dx][i I -u{xY v' {x)dx]d^t{u) 

J E Jr Jr 2 

This motivated the definition given in [6J of a statistical solution of Burgers equation, as a 
family of random processes Ut whose characteristic functions are solutions of ([T]). 

Here we will take for uq a Levy process with negative jumps, that is a cadlag process with 
stationary homogenous increments and negative jumps, and we are really interested only in the law 
of the increments. Hence really defines a probability law on (i?, Ci?)). But u i-^- Jg^u{x)^v' {x)dx 
is not C'{E) measurable, hence this evolving equation for /x makes no sense. One way around this 
is to take a different approach and to work with the Hopf-Cole construction as in [3]. Another 
way was explored in [6^ : it is based on the notion of intrinsic statistical solution, and involves 
working with u{x, t) — u{x, t)dx when b — a gets large : 

Definition 1. An intrinsic statistical solution of Burgers equation is a set {^t)t>f) of probabilities 
on {E,C'{E)) such that for any v e "Dq, 

f I f I , . 1 




dtfit{v) ^ - ^ _hm_^^^y I / {u{x)-- — -/ u{y)dy ] v' (x)dx \ exp[i j uv] dnt{u)) 

In order to work with Levy processes we will need their characterization by means of their 
exponent. 

Definition 2. If u is a Levy process, its characteristic exponent ijj is defined 6i/Va:<y,VAGM, 

(2) ^{exp z[A(u(y) - uix))]} = exp[(2; - x)4>{X)] 

If the process is of finite variance with negative jumps, one can use analytic continuation to 
define its Laplace exponent (f>{X) :— ip{~'iX) as a function on M+. 

In [6] it is shown that if uq is a Levy process of finite variance such that 0o(O) > 0, and if we let 
u evolve according to the non forced Burgers equation, Ut is still a Levy process; its characteristic 
exponent f/'t satisfies the evolution equation dtip = id\{4i'^). 

Let now / be the random forcing, independent on uq. We will denote hy pip its probability 
law. Restricted to events involving a fixed time, /ij^ is the law of a brownian motion defined on M 
by gluing one brownian on and the reflection of another independent brownian on R~ . Such 
a brownian on R is a continuous process with independent homogenous increments null at 0. 

The basic idea now is the following : if ut is a Levy process with probability on C'{E) jit, 
and if we let Ut evolve according to our forced Burgers equation, u{x,t + dt) will be the sum 
of UNF{x,t + dt) and VdtB{x), where unf is the process we would have obtained without any 
forcing, and B(x) is a brownian motion independent of unf- u{x,t + dt) should thus still be a 
process with homogenous independent increments, that is, a Levy process with exponent the sum 

2 

of the exponents. Since the characteristic exponent of a brownian is and since the evolution 
equation for the characteristic exponent in case of no forcing is dt'4^{t, w) — idwip(t, w)'ip{t, w), the 
evolution equation of the characteristic exponent of fi with random forcing should be dtip{t, w) — 
idU^^) + ^. 

One can also get this result in a more rigorous way by working with characteristic functions. 
Let us first define statistical solutions for the forced Burgers equation. Let u(x, t) be a (weak) 
solution of our forced Burgers equation with u{., 0) = uq, u{.,t) e E for t > 0. Let /it denote the 
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law of u{.,t) on C{E). Using Ito's formula, one then gets for each v in T) 
dtfit{v) = / dt{exp[i / u{x)v{x)dx]} dfj.t{u) 

JR 

dt{exp[i / u{x,t)v{x)dx]} dfiQ{uo)diip 
; Jr 

exp [i / u{x,t)v{x)dx)] [ dt[i / u{x,t)v{x)dx] 



fJe 



FJE 

1 
2 



E 



v{x)v{y){f{x, t)f{y, t))tdxdy I dfio{uo)dfiF 
i-u{x)'^v' {x)dx — - I j v{x)v{y)A{x,y)dxdy] exp [i j uv]diii{u) 



Hence our definitions : 

Definition 3. A statistical solution of the forced Burgers equation (*) is a set {iJ-t)t>o of proba- 
bilities on (E,C{E)) such that for any v CzV, 

dtfit{v) ^ 1 ( [ i^u{x)'^v'{x)dx] exp[i j uv]dfit{u) - (l- [ [ v{x)v{y)A{x,y)dxdy] fit{v) 
Je \Jr ^ / JR V^J Jrxr j 

Definition 4. An intrinsic statistical solution of the forced Burgers equation (*) is a set {^t)t>o 
of probabilities on {E,C' {E)) such that for any v S I?0; 

t , . f 



dtfitiv) = ^ -iiT-^+ooy, I / 2 \^^^~b~a ' '^^y^^yi \ exp [i / uv]d^it{u) 

(3) ^ ^v{x)v{y)A{x,y)dxdy^ fit{v) 

One can remark that for A{x, y) = + — jx — yl) and v ^V, jj^^^ v{x)v{y)A{x, y)dxdy = 

Jf^°^ {w{x)^ + {w{—(X)) — w{~x))'^)dx where w{x) = v{y)dy. If moreover v € 2?o, 
!!Ry,Rvi^My)Mx,y)dxdy = ^^w{xfdx. 

A justification for this definition is the proposition, proven in [6], that if u is a homogeneous 
process of finite variance, then 



a— > — oo,b— ^+oo j \ Jjg> 2 \ h — a 



i—u{xYv'{x)dx\e'x:p[i I uv]dfi{u) 



3. Levy solutions 

3.1. Evolution equation. Now we can look for solutions that would be Levy processes with finite 
variance. Such processes are characterized by their exponent ipt, or by their Laplace exponent 0f. 
Let us recall here a useful lemma characterizing Laplace exponents of Levy processes. 

Lemma 1. A function (p : i— > M is the Laplace exponent of a homogeneous Levy process with 
finite variance and negative jumps if and only if (j) is C°° on ]0, +oo[, on 0(0) = and (j)" 
is completely monotonous. 

Moreover, Bernstein's lemma asserts that a fonction g continuous on M+, C°° on ]0, +oo[ is 
completely monotonous if and only if Vn G N, (— l)"^^"-' > 0. Injecting the definition of the 
exponents m in the evolution equation, one gets the following result. 
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Theorem 2. Let {'tjjt)t>a be a family of exponents of Levy processes of finite variance, such that 
VA t I— > iptW is difjerentiable and dti^t is locally hounded. Then 4't characterizes an intrinsic 
statistical solution of our forced Burgers equation if and only if it satisfies the equation : 

w 

Vw e M, dti>{t, w) = idwipit, w)il^{t, w) — 

If the family of exponents has negative jumps, the equation for (f> is 

W 

(4) dt(j){t, w) ^ ~d^^{t, w)4>{t, w) + — 

Proof. Suppose that Ut is a family of Levy processes of exponents (i, w) "0(^7 w) . If u is a fmiction 
in 2?oi the definition of ■0 gives jlt{v) = exp(Jjjj -(/'(i, °° v{y)dy)dx) (as can be seen by considering 
first simple functions). Hence, thanks to local boundedness, dtfj.t{v) = fit{v) /jj dtip{t, v{y)dy)dx 
fitiv) Jjj dtip{t, w{x))dx where we denote as before w{x) = v{y)dy. The right hand side of ([3]) 
is more complicated to get. One needs to compute 

limQ^_oo,b^+oo /e(/r^5("(2^) ~ bhi la Hy)dy)'^v'{x)dx))exp[i J^uvjdfitiu) in terms of i/'- This 
was done in [6] where it is shown that this limit is i{jjg^dw4'(tjw{x))ip{t,w{x))dx)jj,t{v). The other 
term involves 

jj^^^v(x)v[y)A{x,y)dxdy. We have to write it in terms oi w : as noticed before, the specific 
form of A yields JJg^^^v{x)v{y)A{x,y)dxdy = J^v?{x)dx. Finally, we must have for all v in Pq, 
^^dti'{t,w{x))dx = i J^dw'ipit,w{x))ip{t,w{x))dx — ^ ^^w'^{x)dx. This is true if and only if 
Vw e R, dtip{t, w) = id^ipit, w)'ip{t, w)-^. □ 

3 

The stationary solution of equation (|4]) is 4){w) — : this is the exponent of a Levy-stable 

process (which of course is not of finite variance: (j) is not in 0). Speaking informally, this is 
an "invariant measure" of our forced Burgers equation. 

Equation @ can be solved using characteristics. Characteristics for ^ are curves {-^xq (s), ^xq (s)) 
that are solutions of the system 




with initial condition 7^0(0) = xq, Aa;p(0) = (^^{xf)). Each xq corresponds to a different character- 
istics. To find the value of 0(a;, t), one has to find xq, if it exists, such that 7^^^ {t) = x\ if xg exists 
and is unique, then (f)(x,t) = Xj.o{t). 

One can notice that system (I3.1|) describes the motion of a particle in a potential V{x) = 
starting at xq with a velocity (j)o{xo). Thanks to energy conservation, it is thus equivalent to 
A - 7', 73 - 3(7')' = - 30o(xo)^ 7(0) = xo. 

These remarks, and the preceding lemma characterizing Laplace exponents, are the basic in- 
gredients of the next result. 

Theorem 3. Let (j)Q be the exponent of a Levy process of finite variance with negative jumps such 
that 00 (0) > 0. Then the equation ^ with initial condition (/jq admits a unique solution which 
is for all t > the exponent of a homogeneous Levy process with finite variance and negative 
jumps. Such a Levy process is thus an intrinsic statistical solution of the randomly forced Burgers 
equation. 

U 00 = 0, the solution is 4>{x,t) = t^'^ ^ [xt'^Y ~ F^xt'^Y where F is the inverse of x 1— > 
T"). Here V is the Weierstrass function and L02 its half-period. 

We will start the proof by looking at the particular case (f>o = 0. It is instructive to see how 
things go, and it is anyway necessary to deal with it separately. The general case will be covered 
afterwards. 
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3.2. The particular case (f>Q = 0. In this particular case (corresponding to Uq constant), the 
system (13.11) is equivalent to A = 7', — 3(7')^ = a;Q,7(0) = xq. The solution of this equation is 



7(5) = 4:3XqV(uj2 + °i s) where V is the Weierstrass function solution of y — + 1 = (with 

invariants 52 = and 93 = 1), and uj2 = ^J' is the half period of V, satisfying V'{lo2) = 0. 7 is 
1-1 from [0, ^'^^i? 0-12 [ to [xq, +oo[ (see [I]); it is C°° and strictly increasing on ]0, ^'^^ 0^2 [. In order 
to find (j){x,t), one needs to find xq such that 4^a;o7'(ci;2 + t) = x. Hence xq = ^^^^l ^ where 
F is the inverse of a; 1-^ 4:3x^V(uj2 H P is continuous from [0, +oo[ to [0, V443(jj2[ and 

14 ^ 

on 10, +oo[. Near 0, one has A3x^V(u)2^ ^) = x^ + ^ + 0{x^), hence F{x) = i/i- ^ + o(a;^). 

\/34TT ° 

Thus the solution of jH with initial condition 0o = is (j){x,t) = t^^^{xt'^Y — F{xt'^)^; the 
development of F near ensures that x 1— > 0(a;,t) is on [0, +oo[. The proof that 0" is 
completely monotone will be done for the general case. 

3.3. General case. As a preliminary remark, since is completely monotone, 0o > on K+, 
hence (j>Q is nonnegative and increasing on R+. Moreover, if there is an interval I — [0, a[ such that 
4>q{w) = OVw G J, then 0q = on /. But is decreasing (because 0q^^ < 0) hence 0q = on M+. 
Therefore 0o = on M+ or 0o > on ]0, +oo[. The particular case (f)Q — Q has been covered in 
the preceding section (except for complete monotonocity), hence we will assume here that (/)o > 
on R+*. 

Everywhere in this proof 0(a;, t)' ,(f>{x, t)" , (j){x, t)^"' will refer to partial derivatives with respect 
to x. 

As noticed before, the solution 'jxoit) of the system p.ip is the trajectory of a particle in 
a potential V{x) — — starting at a^o with a velocity (I)o{xq). Hence jxo{t) is a continuous 
increasing function of xq. More precisely, if a;o = 70 (0 = 0, else if xq ^ 0, 



7x0 (t) = x^t= j 

J XQ 



dy 



y/y^ -xl + 3c^l{xo) 

The function h : {(x,xq)/0 < xq < x} R+, defined by h(x,xo) = , is C°° 

JXo ^j,3_3.3_|_302(j.^-| 

{(j)Q{x{)) > 0, therefore this integral is well defined); xq i— > h(x,xo) is a diffeomorphism from ]0,x] 
into hence for alH > jxo{t) = x defines uniquely xo as a C°° function of {x,t) on ]0, +oo[^, 
satisfying limj;^o xq{x, t) — 0. 

Hence the equation (g]) admits as a solution for all t, (j){x,t) = y/ x^ — xo{x, t)^ + 30§(xo), 
(j){0, t) = 0. X ^ (j){x, t) is C°° on ]0, +00 [. We still have to check that (j)" is completely monotone 
and that is at 0. 

Let us prove that 0" is completely monotone. This proof will also work for cj)Q = 0. The keypoint 

is to notice that ^((/.(7,„(t), i)) = A,„(t)0'(7,„ (i), t) + ^(7,^ (i), i) = 0</''(7xo + ^(7.0 i)- 
Hence by deriving twice equation (jH) with respect to w, and by setting w — 7^^ (t) one gets 

|(0"(7.„(t),t))+30V'(7.„(i),i) = l 

Thus, for every xq, 1 1-^ 0"(7xo (Oi ^)) satisfies a differential inequality of the form ^ + 3/it > 
(where / = 0'(7a;o(i),t)). Since 0q > 0, we deduce 0"(7a;o(t), i) > for all t > 0. Since it is true 
for aU Xo, 4i"{x, t) > for all a; > 0, t > 0. 

One can now proceed by induction. Suppose (— l)"!/)*^"-* (x, t) > for all 2 < n < N. Then by 
deriving ^ N + \ times one gets 

N / \ 

^''^'HlxMht))HN+2)f^^''+'\j.„{t),t))=-W 



Now the induction hypothesis guarantees that the right hand side of this equation is of the sign 
of (-1)^+1. Hence {-l)'^+'^(j)^^+^\j^„(t),t) verifies ^ + {N + 2)fu > 0. One concludes as before 
that i-l)^+^4>^^+^\x,t) > 0. 
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Hence (f){x, t) is completely monotone. 

Let us now prove that is at .x = : wc need to know how xq behaves near x = 0. In 
order to do this we will consider the solution of — 3y'^ + C{xo) = 0, y{0) = xq, where C{xo) = 
3(^o(a;o)^ — Xq- Its explicit expression depends on the sign of C. We know that lima;^oa;o = 0. 
Hence there arc really two cases that have to be dealt with : (l>'o{0) > and ^'q{0) = 0. 
. //</.[,(0) = 

(f)o is at 0, hence (/>o(a;o) = k'^x^ + o(xq) where fc is a constant, and C{xo) = 
—Xq + 3k^XQ + o{xq) is negative when xq is small enough. Hence 7xo(i) = AV{bt + c) 

where = -4C(xo), b = V is as before the Weierstrass function with invariants 

92 = 0,93 = 1, and c verifies AV{c) = xq. 

Prom the expansion V{w2 + e) = 4~3 + 3.4~3e2 _|. o(e^) we deduce 
A = 43xo(l - fc-^xo) + o(x^) 

b = 




/ 43 „ 

W2 + \l —k + o(v^) 

Ixoit) = xo + {kH+ —)xl + o{xl) 
XQ{x,t) = X- [h^f + — ).T^ + o(.T^) 



= ^--^^^^±^^^ = {k^ + \)x-+o{x^) 

Hence is at x = 0. 
• If<i>'Q{Q) = a>Q 

Now C(xo) = Zo^Xq + kxQ + o(xq), where fc is a constant, is positive when xq is small 

A 



enough. Hence 7^0 (^) = AQ{bt + c) where A^ = 4C(.to), b = y^, and Q is now the 
Weierstrass function with invariants g2 = 0,93 = —1; and c must verify AQ{c) = xq. 
Q is related to V by Q{z) = —'P{iz) for all in C. Q admits a real zero z\, with Q'{z\) = 1 
and Q"{z\) = Q"'{zi) = 0. It is not difficult to get the expansion Q{zi + e) = e+^+o{e'^), 
and : 

A = (12a2)ix|(l + ^xo + o(a;o)) 
,axo^l/, k , 

/ Xq ,1 _8,12^3 12^3 4 4 



<^{x,t) = y ^^ = ^^Vl + ^2ar + o(a;) where = i2a^(i+atV 

Hence ^ is here also at x = 0. 
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